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SUMMARY 


This paper deals with theoretical analysis of flow induced instability 
force due to labyrinth seal. That is, a approximate solution is given for the 
partial differential equation representing the flow in labyrinth seal and it is 
compared with the finite difference method which was proposed in the previous 
report in order to verify the accuracy of both methods. Then the effects of 
difference of inlet and outlet pressures of the seal, deflection of pressure and 
mass flow from the steady state, rotor diameter, seal clearance, seal interval 
and seal number on the flow induced force of the seal are investigated and it is 
known that some of these factors are very influential on the flow induced force. 


INTRODUCTION 


In the previous report (l), the fxmdamental equation of flow in the labyrinth 
seal was derived by considering the effect of the variation of gland cross section. 
The equation is numerically solved by using the finite difference method. Then 
the spring and damping coefficients of the labyrinth seal are calc\ilated, and by 
using the result, the stability of a rotor system having labyrinth seal is 
discussed on this coefficient by using the energy concept. Furthermore, experi- 
ments are executed to observe the flow pattern in the gland and to study the 
characteristics of the flow induced forces in the labyrinth seal. 

However, further theoretical investigations are required for the accirracy of 
niimerical calculation and more detail description of labyrinth seal behaviour. 

Then in this paper the equation is solved by another method and the calculated 
results are compared with the solution of finite difference method. After that 
influences of deflection of steady state, pressure difference between inlet and 
outlet, labyrinth seal radius, seal clearance, seal pitch, number of seal cham- 
ber, seal strip height and divergence and convergence seal on the flow induced 
force, phase angle and leakage flow rate are investigated. 


SYMBOLS 


t time 

w = Rsf 

Rs radius of labyrinth seal 

p pressirre in seal chamber 
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c peripheral velocity in seal chamber 

f density of gas 

q mass flow rate in axial direction 

y angle from x axis 

f cross section of seal chamber 

length of strip pitch 
h strip height 

§ radial clearance of seal 

X' friction shear stress of stator surface 

I friction shear stress of rotor surface 

U' length of acting surface of shear stress 

U" length of acting surface of shear stress 

R gas constant 

T absolute temperature of gas in seal 

7' specific weight of gas 

flow coefficient 
n ratio of specific heat 

CO rotating speed of rotor 

60 n critical speed of rotor system 

u=(ORs peripherical velocity of rotor 
K coefficient of viscosity 

U coefficient of kinematic viscosity 

Re Reynold's n\miber 

}\ friction coefficient (stator) 

X* friction coefficient (rotor) 

subscript * steady state 


( stator) 
(rotor) 


206 



EQUATION OF MOTION 


Derivation of Equation of Motion 

Derivation of the equation which describes the labyrinth seal flow is 
followed the previous report (l), where the following assumptions are taken and 
schematic flow is shown in Fig.l; 

(1) Fluid in labyrinth seal is assumed to be gas and its behavior is assumed to 
be ideal. 

(2) Temperature of the fluid in the labyrinth seal is assmed to be constant. 

( 3 ) Cross section area of the seal gland is assumed to be constant in spite of 
the deflection of rotor and time derivative of cross section area is only consi- 
dered. 


(H) Change of flow state in the gland is assumed to be isentropic change. 

Thus the fundamental equations with respect to continuity and momentvun for the 
ith seal element are as follows; 


Q(Pih) 

9t 



9CfiU) 

dVT 


i ( — 0 


( 1 ) 


dvr 


4 




( 2 ) 


As temperature of 
is given as; 


each gland is constant, 

?t- nRj 


the equation of state in the ith gland 


C3) 


where 


ii= A f 


Next, denoting the gas flow velocity through labyrinth nozzle (seal clearance) 
S, the relation between axial flow rate and pressure in the gland is given by 
the thermodynamical energy equilibrivim condition as. 


S = 



K1 

n- 


Pi^i 



(4) 


where n is isentropic exponent, v specific voliune and subscripts I and 0 
mean inlet and outlet of nozzle, respectively. Assuming that the flow is adiabatic 
change, i.e. Vo =(Pi / Po)^xVj, the flow rate through the nozzle G [kgf/sec] is 
given as; 

.. iC.S .F-ii/ ro t: Ir,^ H KII/ I 

(5) 


I r* • ^ 3 — — 
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where F is nozzle cross area (area of seal clearance) and subscripts (i,0) 
correspond to the gland number of the labyrinth seal. Applying the above rela- 
tion to flow in the labyrinth seal, we assume that the labyrinth seal is a 
series of the nozzle, pressure difference between each gland is small and enthal- 
py of each gland is nearly equal. Then leaJtage flow rate through labyrinth seal 
G (kgf/sec) is given by using the relation 
P, V, =P, V, = • ■ • =p, V. as ; 



( 6 ) 


Applying Eq.(6) to the flow through the ith seal element, mass flow rate for 
unit area is given as; 


_5l 

hi F-J 


0^ 




Rewriting Eq.(T), the following equation is obtained; 

£ 


(7) 


where 




A" 


_ »RT 


( 8 ) 

(9) 


Linearization of Fundamental Equation 

For the linearization of the equation, the perturbations of 'pressure, 
peripheral velocity and flow rate from those of the steady state and also mean 
value and perturbation of seal clearance are introduced as 

+ C^=C*LCl + t.J 

where g,^,c^^and q„^are pressure, peripheral velocity and axial flow rate of 
steady state in the ith gland and , *^j^and l^^are nondimens ional pertiirbation 
terms of pressure, peripheral velocity and axial flow rate. Also and hr 
are mean value and perturbation of seal clearance. 

Assuming that the rotor is whirling along elliptical orbit, is repre- 
sented as 


= coswt cos<f - ^SINwtSlNf 

0*1 


(11) 
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Then the cross section area in the ith gland is represented as 

As the change of state in the gland is isentropic, the following relations are 
obtained; 

-fe = -^if = A (13) 


fz~ n 


therefore 


_9H.= n-S.^ 

9t A 9t 

Substituting Eqs.(lO) - (lU) into Eqs.(l) and (2), and denoting the time 
derivative and spatial derivative {• ) and respectively, the 

following equations are obtained for the ith seal element ; 


T*u ► , 

n HRs ' 

n n-r-r»2 tjX 


Cif'fiii, f,/ , 




i,V.i 2Ri,&. 

= (Rl-p,t„)[ai.,cosa)tCOSf+bi.,SmaASin^f j 

b^i. ( R,ir pyj{ai,coswtcosf + bi,si nwts'mf ) 

-f y [-aLtDSimutcoSf + bicucoSLutimf-aLtiii^Sina/tSinf +bbt|U)Cosa^tsinf J 

( 

fe I X _l_ 1 ^ RT . 0^< i/ ^ 4- [.2 ^ R L'b / '1^ V 

,[ ^RTAV^u+.p*t , ^RT/i^SLUl , Xlur ^ , 2 ]^ 

+ XiC.vU(+XlU"(U-6,jj<t,- ti-r h; 

+ XlOtUi. >^iUl f , p ^2 , ^RTA^/r-2 -n rl pJr^i _» Ti^.. _9 _ . 

- i^^^-JiLF,J(aiCOSwtCoS<f+USinu>tSinfl-hj-(-a-t^Siniutcc:S‘f 

-f biiJ^C^SivtSii^i’-ai-t-iWSina^CoSfi-bitilJUCoSwtSincf) , 
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For the n stage labyrinth seal, Eqs.(l5) and (l6) are represented in matrix form 

Till + V Ul t Aui = +6ut )+ CUcoS(<f-ajt)-hiiCoSi^+a)t) 

+ bjtCoSC^f-o;!^ + G2[-^USin(tf 

+ b*5in(f -ort^ + G^[Qoi(^S(f+lOt) + 6a(^K^-iui)--0hCjOS(}{^ + tx;t) 

+ 0b coS(if _ s in (4> + a;t) 1 9aSw C^-ivt) + 0bSin + o;t; 

+ 6bS in (cp_£tjt}j (17) 

where T-^o= ?o (inlet pressure) and Fun = Pn (outlet pressure). 


ANALYSIS 


Analysis of Steady State 

Setting the perturbation terms ^j_and 'll zero in Eqs.(l5) and (l6) in order 
to obtain the pressure and peripheral velocity of each seal stage in steady 
state , 

«i+l' ^1 ' ^ <18) 


q-_L-iOU-- <i.(V . + t'u.' - x"u." = 0 

^1+1^ 1 1-1 1 1 


( 19 ) 


As the multi labyrinth stage is considered as a serial nozzle, the axial flow 
rate in the steady state is obtained as follows 



(20) 


where P'l-i is the pressure of (i-l)th chamber which is taken into account of 
draughty flow. The pressure is given by Ref. (3) as; 


B:=?L.+(po-p«)(i-A05ptp 


u _ 2 Aj, 

Alt I 


( 21 ) 


where A^ is determined by the shape of labyrinth seal; 



^Cc+tan8 


( 22 ) 


where Cc, is coefficient of vena contracta. 

From these equations, state variables i.e. the pressure and flow rate satisfied 
the equilibrium condition, is obtained by iterating Eqs.(l8) and (l9). The 
peripheral velocity is obtained by using Eq.(l9) in the same manner. 
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Analysis of Dynamic Force 


Referring to the left term in Eq..(l7), the solutions are assumed as follows 


H = ?^i + Rci { A+L5inC4>+a;t)+B+tC0S(f 

Cl= Citi+C^L[EnSin (tf+ FHC^S(tp+a^t;tE.i,Sin (f-a/t; t 

Representing the above coefficients in matrix form as 


X V 

Ati 


fA-A 


B+r 


fB-.] 


, 0^2= 

E-. 

. fr 

F*i 

11 

F- 

A+n 


A-. 


B+n 



Etn 


E-n 


Fi-n 


F-n 

0 


J 


Lo J 


.0 J 


the solution is written as; 


(23) 

(24) 


U = diS]r)(^+ajt)+fitC0S(^^-UJt)-\-cJl2Sir\C^-u)t)+^2COS(ir~iVt) (25) 


Substituting Eq.(25) into Eq.(lT)j and separating it to the elements cos(^ + a>t), 
sin(9’+^'t), cosC^’-tOt) and sin(y-cot), the coefficients are represented as 
follows ; 

(jOlTci^i + y/ojti + A^i = (Bi(X* t +G fb* + G?0b 


- (jJFA "\^i t Aq^i + (I746A+ Ga b» + 

- (t)To^2 +\/q^2 t A/^2 = + HI 3 Oai~ H?b* + Hl75b 

coTi^2 -yj^z t AQ^2=*-H}2(Z^f + H1404+ Hk bii + Hl8 0*> 


( 26 ) 

(27) 

( 28 ) 
(29) 


where G, /sy Ga and iHi ~ iHg are coefficient matrices which is composed of the 
matrices [g » £, > |p and(Q . Solving Eqs . ( 26 ) , ( 27 ) and Eqs.(28), (29) simultaneous- 

Cjil . i^l ’0^2 and Pz are obtained. Then velocity and pressure in the gland is 
calculated by substituting the coefficients into Eq.(23), and the flow induced 
force for whole seal is obtained by integrating the pressure for each gland. 
Representing the x and y direction force by Px.and , 


Px=-RsiffR,^JcosWV 


The flow induced force F and phase angle from x axis are described as 




<p = tan-' -B 


(31) 


Spring coefficients and damping coefficients are obtained from Eqs. (30) and (3l) 
in the following matrix form; 


211 



( 32 ) 



0 

X ' 


^ V 

CxxCxU Q 


t' • N 

X 

KptKjij 

M 

+ 

Cjx 


9 

0 

Kexx Kexj 

0x 


Q Cexx Cpxa 
Cejx Ceaa 


Bx. 

KeyxKei^a 

/ 

1 ^ 





When Eqs.(30) and (3l) are rewritten to Eq. (32) in the above, the following form 
is used as , 

^iiCl*C0Sa)t-t- ^sihuCoSivt ^^^Q^cosuJt 

■^^^iBaSincvt + ^^^^hCoSwt + ^^^OhSincot 

Using the above expression, the following relations are obtained between coeffi- 
cients; 

Kxx.== Kflxx= Kesa Kxd=^-K3x ^ 

Cxx=Casj Cexx= Cgij^ Cxn = ~Cjx Csx!j="(^6iix 

NUMERICAL EXAMPLE AND DISCUSSIONS 


The labyrinth seal having three teeth, which is shown in Fig. 2, is used as 
a numerical model. Details of the labyrinth seal is shown in Table 1 and for 
the parameter survey each element of seal dimension is independently changed in 
order to investigate the influence of its seal size on the flow induced force. 

In order to evaluate the seal force, circular orbit (a^^ = b^ in Eq. (IT)) are 
used as a whirling of labyrinth seal, and its radius a^ =0.1mm. For the calcula- 
tion the following equation is used as the friction coefficient between rotor 
and stator; 
for stator surface; 


D ' - D( ~ 

Kcu 

where D is wetted perimeter length considered the rotating Reynold's Number 

Ret < I 200 , \'i=U/Rei , Rei 1 1200 , o. 3/^4- 


for rotor surface; 

- Ptu-w 


f 


Ret < 1200 , X'i = 1200, At = 0. 


Two methods are used for the numerical calcination, one is the method which is 
shown in previous section, the other is one which is the finite diference method 
that was proposed in Ref.(l). At first calculated results for these two methods 
are compared with and as the result it is known that these two results agree. 
Then parameter surveys of labyrinth seal are carried out to see the influence 
on the flow induced force which exits the instability. 

In this paper critical speed of a rotor system is assumed to be 5000 rpm, 
and rotating speed of the rotor is represented in nondimens ional expression as 
^/Ct)n(<j0 (rpm) : rotating speed of the rotor). F and FAI shown in figures mean the 
absolute value of the force induced by labyrinth seal and its phase angle from x 
axis . 
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a) Accuracy of Two Calculation Method 

Fig. 3 shows the spring coefficients and damping coefficients which are 
obtained by two methods i. e. approximate method and finite difference method. In 
the figure, full line and broken line show the coefficients for inlet pressure Pq 
= 196.13(KPa) and ll+T.l(KPa) respectively and outlet pressure Pn=98. OT(Kpa) which 
are obtained by the approximate method and doted line shows the coefficient for 
inlet pressure Pq = 196 . 13(KPa) and outlet pressure P^=98.0(KPa) which is obtained 
by the finite difference method. Also sign (-) on line means negative value. 

From this figure full line and dotted line coincide very well so each calculation 
method has good accioracy. 

b) Influence of Deflection of Steady 

In this paper pressure, peripheral velocity and flow rate in steady state 
are iteratively calculated, however, Kostyuk and other authors obtain them by 
another method and there are some discrepancies between both steady state 
values. Moreover, steady state values in field and the calculated values may 
have discrepancies. So it is important to investigate the influence of deflec- 
tion from steady state on the accuracy of the induced force. Then Fig.U shows 
the flow induced force and its phase angle when pressure is a little deflected 
from steady state where the number on lines correspond to the n\miber in Table 3. 

The calculated data are shown in table 2. From the calculated result, it is 
known that even if the flow rate is deflected 20^ from the steady state, flow 
induced force and phase angle are not so influenced, and these curves become like 
one line, so this figure is neglected. On the other hand if the pressure is 
deflected 0.5^ from steady state value, they are influenced very much. Also it 
is known that the induced force is not influenced by rotating speed of rotor but 
if it is separated to component, influence of rotating speed is remarkable as 
shown in Fig. 5. From the above discussions it is concluded that in the calculation 
of flow induced force, steady state values must be calculated with high accuracy. 

c) Influence of Pressure Difference Between Inlet and Outlet 

Fig. 6 shows the influence of pressure difference between inlet and outlet 
on the flow induced force when inlet (entrance) pressure is kept constant value 
(209 atm) and outlet (exit) pressure is changed. The figure shows that as the 
outlet pressure becomes large, i.e. the pressure difference between inlet and 
outlet becomes small, the flow induced force becomes small. As the flow rate 
decreases, the induced force proportionally decreases. 

d) Influence of Labyrinth Seal Radius 

Fig. 7 shows the influence of labyrinth seal radius on the flow induced 
force, where the seal radius is valiable and other dimensions are fixed at 
constant values. Here flow rate means total flow rate as shown in the expression; 


From the figure the flow induced force proportionaly increases as the seal 
radius increases. Therefore, the flow induced force due to labyrinth seal can be 
easily evaluated by the similarity law, when the seal radius is known. 

e) Influence of Seal Clearance 

Fig. 8 shows the relation between seal clearance and the flow induced force 
and the leakage flow rate (leakage mass flow rate per unit length), where as the 
clearance increases, the induced force F decreases and the flow rate increases 
proportionally. 
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f) Influence of Seal Pitch 

Fig. 9 shows the influence of seal pitch on the flow induced force, when 
strip height is kept constant. From the figure it is known that as the seal 
pitch increases, the induced force increases, but the leakage flow rate q does 
not so much changed. Thus the large pitch is not good from the view point of 
the flow induced force. 


g) Influence of Number of Seal Chamber 

Fig. 10 shows the influence of chamber number on the induced force and the 
leakage flow rate, where inlet and outlet pressure is constant and each chamber 
dimension is also same. From the figure the induced force increases like hyper- 
bolic as the chamber number increasing. On the other hand leakage flow rate 
decreased like exponentially. Therefore optimum number of chambers may be 
determined when inlet and outlet pressure are known. 


h) Influence of Seal Strip Height 

Fig.ll(a)(b) shows the influence of seal strip height. In this case the 
induced force, phase angle and leahage flow rate do not influenced so much. But 
if this result is represented in the spring and damping coefficients, the influ- 
ence due to the strip height is remarkable. Fig ll(c) shows the spring coefficients 

k and k 
XX xy 


i) Influence of Divergence and Convergence Type Labyrinth Seal 

Fig. 12 shows a comparison of seal types, that Is, comparison of straight 
through type, convergence type and divergence type. In this calculation there 
is little discrepancies of the flow rate q, so its figure is neglected. From 
Fig. 12(a), it is known that the flow induced force for the convergence type is 
the smallest of all and divergence type is the largest inspite of the same 
leakage flow rate. From Fig. 12(c), the flow induced coupling stiffness k:tx is 
similar characteristics, but convergence type is excellent from the viewpoint 
of stability. 


CONCLUSION 


In this report the flow induced force due to labyrinth seal which is some- 
times the cause of instability is studied for the special model. That is, the 
labyrinth seal of shroud of steam turbine is used as an example and its flow 
induced force and stiffness and damping coefficients are theoretically calcula- 
ted by two method. Then influences of deflection of steady state, pressure 
difference between inlet and outlet, labyrinth seal radius, seal clearance, seal 
pitch, number of seal chamber and seal strip height on the flow induced force 
and leakage flow rate are investigated. As the results, it is known that the 
above factors are very influential on the flow induced force. However further 
theoretical and experimental investigation are required in order to obtain a 
more precise description of labyrinth seal characteristics. 
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TABLE 1 


LABYRINTH SEAL DIMENSIONS FOR NUMERICAL EXAMPLE 


SEAL RADIUS 

Rs 

1 

(m) 

ROTOR NATURAL FREQUENCY 

UJn 

5000 

(rpm) 

STRIP HEIGHT 

h 

10 

(mm) 

SEAL CLEARANCE 

s 

1.3 

(mm) 

STRIP PITCH 

1 

30 

(mm) 

NUMBER OF CHAMBER 

N 

2 

SEAL TYPE 


straight 

-through 



type 


TABLE 2, - GAS CONDITION FOR NUMERICAL EXAMPLE 


OPERATING LIQUID 

air 

LIQUID TEMPERATURE 

T 

784.15 

(K) 

ENTRANCE PRESSURE 

Po 

20482.0 

(kPa) 

EXIT PRESSURE 

Pn 

19678.4 

(kPa) 


TABLE 3. - DEFLECTION OF PRESSURE AND FLDW RATE FROM STEADY STATE 



PRESSURE OF EACH STAGE 

(kPa) 

FLOW RATE 
(kg s/m*) 

PO 

PI 

P2 

P3 

STEADY STATE VALUE 

(1) 

20482.0 

20166.3 

19936.7 

19678.4 

6.8443 

DEFLECTION OF 


20482.0 

20267.2 

20036.4 

19678.4 

6.8443 

PRESSURE P1.P2 

BB 

20482.0 

20065.5 

19836.9 

19678.4 

6.8443 

DEFLECTION OF 


20482.0 

20166.3 

19936.7 

19678.4 

7.5287 

FLOW RATE 


20482.0 

20166.3 

19936.7 

19678.4 

6.1599 
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Figure 3. - Accuracy of two calculating methods and influence of pressure 
difference between inlet and outlet. 
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Figure 4. - Influence of pressure deflection from steady state on 
the flow induced force. 
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Figure 5. - Influence of pressure deflection from steady state on 
the cross coupling stiffness. 
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